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Abstract 
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1 Introduction 



Perturbative QCD predictions [l] are important to test QCD and to measure its funda- 
mental parameters (e.g. the strong coupling as), to estimate Standard Model backgrounds 
for new-physics signals, and to understand the interplay between perturbative and non- 
perturbative phenomena. 

Calculations at the leading order (LO) in perturbation theory are straightforward: the 
matrix elements at tree level can be easily computed by using helicity techniques and colour- 
subamplitude decompositions and can then be numerically or analytically integrated. 

Nowadays, also next-to-leading order (NLO) calculations are feasible in a direct man- 
ner, as witnessed by the accelerated production rate of new NLO computations 0. This 
achievement is the result of much theoretical progress in the last few years. This progress 
regards both efficient techniques for the evaluation of one-loop matrix elements |4| and com- 
pletely general algorithms f||-§] to handle and cancel infrared singularities when combining 
tree-level and one-loop contributions in the evaluation of physical quantities. 

In contrast, the calculation of jet observables at next-to-next-to-leading order (NNLO), 
although feasible in principle, remains a challenge in practice. Indeed, no computations 
have appeared that involve more than two kinematic variables @. The general obstacles 
that have to be overcome to perform NNLO calculations are the evaluation of two-loop 
amplitudes and the cancellation of infrared singularities between real and virtual contri- 
butions. There is only one QCD amplitude that is known at two-loop order, namely the 



electromagnetic form factor of the quark [10, 11]. The extension to higher loops of cutting 



techniques 0] has recently led to the computation of the two-loop four-gluon amplitude in 
the simplified case of N = 4 supersymmetric QCD [0. Other approaches for efficiently 



calculating higher- loop amplitudes taken by various authors include string theory |T3| , first 



quantized [14]] and recursive approaches. Regarding the cancellation of infrared sin- 



gularities, much work still has to be done |T6| to extend the NLO algorithms to higher 
orders. 

In the case of calculations at NLO, an important ingredient of the theoretical progress 
mentioned above has been the complete understanding of the factorization properties of 
tree-level JTJ and one-loop JT7], [5], ^ amplitudes in the soft and collinear limits. An analogous 
understanding at NNLO is important as well. In the case of tree-level amplitudes, the 
factorization structure in the double-soft |18| and double-collinear 0] limits is known. 



The collinear limit of one-loop amplitudes [^] is also known. The singular behaviour of 
two-loop amplitudes has not been investigated so far. 

In this paper we discuss the structure of infrared singularities in on-shell QCD ampli- 
tudes at two-loop order and, in particular, we present a universal factorization formula that 
gives the coefficients of the e-poles of the dimensionally regularized divergences. Our results 



have been obtained by improving and extending the coherent-state approach of Ref. |2T[ to 
higher infrared accuracy. More details on the derivation of the results presented here and 
on the singular behaviour of all-loop amplitudes will be given elsewhere. 

The outline of the paper is as follows. In Sect. El we introduce our notation and, in 
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particular, we point out the subtleties related to the regularization scheme of the infrared 
divergences in the general context of dimensional regularization. Section |^ is devoted to a 
review of the known results on the singular behaviour of one-loop amplitudes. Our general 
factorization formula for two- loop amplitudes is presented and discussed in Sect. f| and 
applications to some particular cases are considered in Sects. |5] and ||. At present, we are 
unable to explicitly compute the coefficient of the single poles in general terms, but we can 



do that in the particular case of amplitudes with two quarks, as shown in Sect. |5T2| . Our 
results are summarized in Sect. |7|. 



2 Notation 



2.1 Dimensional regularization 

In the evaluation of loop amplitudes one encounters ultraviolet and infrared singularities 
that have to be properly regularized. The most efficient method to simultaneously regu- 
larize both singularities in gauge theories is to use dimensional regularization [p2}-|2"4j|. 



The key ingredient of dimensional regularization is the analytic continuation of loop 
momenta to d = 4 — 2e space-time dimensions. Having done this, one is left with some 
freedom regarding the dimensionality of the momenta of the external particles as well as 
the number of polarizations of both external and internal particles. This leads to different 
regularization schemes (RS) within the dimensional-regularization prescription. 

The RS known as conventional dimensional regularization (CDR) is, in a sense, the 
most natural scheme because no distinction is made between particles in the loops and 
external particles. All particle momenta are <i-dimensional and one considers d — 2 helicity 
states for gluons and 2 helicity states for massless quarks. Nonetheless, other RS f\ are 
known to be useful for field theoretical reasons (e.g. to explicitly preserve supersymmetric 
Ward identities [25]) and/or because of practical simplifications in actual calculations M. 



This brief overview of dimensional regularization is important to discuss the singular 
behaviour of loop amplitudes, because only the leading singularities are RS-independent. 
In general, the coefficients of all the other singular terms as well as the finite parts do 
depend on the RS. The RS dependence regards both ultraviolet and infrared singularities, 
but it affects them in different ways. 

In principle the regularized ultraviolet divergences have to be removed from off-shell 
Green functions via renormalization of wave functions and coupling constant. This leads 
to the introduction of the running coupling as(n 2 ), whose definition in terms of the bare 
coupling depends on the regularization (and renormalization) scheme. The ensuing scheme 
dependence in renormalized on-shell amplitudes is nonetheless trivial. It simply amounts 
to (and can be controlled by) an overall perturbative shift in as (cf. Sect. [2.2|) . This 
perturbative shift has no effect on tree-amplitudes and is universal in loop amplitudes, 



'''The RS that are mostly used in one-loop computations are the 't Hooft and Veltman scheme ]2^ j, the 
dimensional-reduction scheme pj| and the four-dimensional helicity scheme pfl . 
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that is, independent of the number and type of particles. 



Infrared singularities behave differently with respect to dimensional-regularization pre- 
scriptions. Although one can keep track of their RS dependence in general terms [27|, to 



relate loop amplitudes in different RS one has to introduce proper transition factors that 
depend on the number and type of external partons. In the context of one-loop amplitudes, 
this feature was pointed out in Ref. [p8|| , where the (one-loop) transition factors between 



the main schemes in current usage were explicitly computed. Beyond one-loop order, care 
has to be taken in the operational definition of some RS. We do not explicitly consider this 
problem, but in Sects. |]-|6| we discuss the RS dependence of two-loop amplitudes within 
the general class of schemes specified at the end of Sect. O. 



Note also that, in order to fulfil unitarity in calculations of physical quantities, the RS 
dependence of loop amplitudes has to be consistently matched to that of tree amplitudes. 



This issue is discussed on quite a general basis in Ref. |27| and is no longer considered in 
this paper. 



2.2 Renormalized amplitudes 

We consider amplitudes M. m that involve m external QCD partons (gluons and massless 
quarks) with momenta pi, . . . ,p m and an arbitrary number and type of particles with no 
colour (photons, leptons, vector bosons, ...). Note that, by definition, we always consider 
incoming and outgoing parton momenta in the physical region, i.e. any pi is massless and 
with positive-definite energy (in particular, p^ ■ pj > 0). The amplitudes are denoted by 
•M-mivi, ■ ■ ■ iVm) ( or 5 shortly, M. m ({p})) and the dependence on the momenta and quantum 
numbers of non-QCD particles is always understood. Besides that, the unrenormalized 
amplitude M. m depends on powers of a$ (Mq, where «g is the bare coupling constant and 
Ho is the dimensional-regularization scale. 

We find it convenient to use renormalized amplitudes. These are obtained from the 
unrenormalized ones by just expressing the bare coupling a$ in terms of the running 
coupling as(/i 2 ) evaluated at the arbitrary renormalization scale /i 2 . In order to avoid 
renormalization-scheme ambiguities, we always consider the running coupling as(fi 2 ) as 
that defined in the standard MS scheme. Thus, we use the following expression 



$ S e = a s (/i 2 ) /i 2e 



«s(^ 2 )^ + «i(/i 2 ) +0(altf)) 



e ^ ' \e 2 2e, 
where /3q, (3\ are the first two coefficients of the QCD beta function 

„ UC A -AT R N, nCl-WC A T R N f -6C F T R N, 

A = vT, ' ft = w ' (2) 

and S e is the typical phase-space volume factor in d = 4 — 2e dimensions (je = — ^(1) = 
0.5772 ... is the Euler number) 

S e = exp[ e (ln4vr + ^(l))] . (3) 



Note that, as recalled in Sect. |2.1|, unrenormalized amplitudes contain some RS depen- 



dence of ultraviolet origin. Using the MS renormalization scheme, we have eliminated this 
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harmless dependence from renormalized amplitudes. In practical terms, our renormaliza- 
tion procedure works as follows. If the bare amplitude is evaluated in CDR, one simply 
replaces according to Eq. (H). If a different RS is used, the minimally-subtracted sub- 
stitution in Eq. (|1|) still applies in terms of an RS-dependent coupling ag a (/z 2 ). To obtain 
the MS amplitudes, one then has to perform the additional substitution: 



a 



RS, 



1 + di a s + d 2 «s + C( a s 
where the coefficients di, d 2 depend on the RS and are finite for e — > 0. 

The MS-renormalized amplitude has the following perturbative expansion: 

/ a s (/i 2 )\ 



(4) 



2tt 



+ 



2vr 



(5) 



where the overall power q is half-integer (q = 0, 1/2, 1, 3/2, . . .), in general. This equation 
fixes the normalization of the tree-level (.M^?), one-loop and two-loop (Adff) co ~ 

efficient amplitudes that we use in the rest of the papeiQ. These coefficient amplitudes are 
infrared-singular when e — > and behave as 



-Ml n) 



2d 



+ 



(6) 



where the dots stand for e-poles of lower order. The coefficient of the poles are still RS- 
dependent. The results presented in the following sections are valid in CDR and in any 
other RS that is consistently defined as ultraviolet regulator. By that, we mean in any other 
scheme whose minimally-subtracted coupling is related to the MS coupling via Eq. (f|) in a 
universal way (i.e. with coefficients di, d 2 that do not depend on the computed amplitude). 



2.3 Colour space 

The colour structure of QCD amplitudes can be handled in two different and equivalent 
ways. One method consists in projecting out the amplitude onto a particular set of non- 
orthogonal colour vectors, thus extracting colour subamplitudes ||. The other method 
works directly in colour space. We prefer to use the colour-space formalism and, in partic- 
ular, the same notation as in Ref. ||. 

The colour indices of the m partons in the amplitude Ai m are generically denoted by 
Ci, . . . , c m : Cj = {a} = 1, . . . , — 1 for gluons and Cj = {a} — 1, . . . , N c for quarks and 
antiquarks. We formally introduce an orthogonal basis of unit vectors {\c\, . . . , c m )} in the 
m-parton colour space, in such a way that the colour amplitude can be written as follows: 

M%-' Cm {pi,...,p m ) = (ci,...,C n \M m {Pl,---,Pm)) ■ (7) 

■^Precisely speaking, Ad^m is not necessarily a tree amplitude, but rather the lowest-order amplitude for 
that given process; Mm ,-M-m are the corresponding one-loop and two-loop corrections. For instance, in 
the case of gg — > 77, A4^ involves a quark loop. 
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Thus \JA m (pi, . . . ,p m )) is an abstract vector in colour space and, in particular, the square 
amplitude summed over colours is: 

\M m ({p})\ 2 = (M m ({p}) | M m ({p})) . (8) 

Colour interactions at the QCD vertices are represented by associating a colour charge 
Tj with the emission of a gluon from each parton i. The colour charge Tj = {T"} is a 
vector with respect to the colour indices a of the emitted gluon and an SU(N C ) matrix 
with respect to the colour indices of the parton i. More precisely, its action onto the colour 
space is defined by 

(ci, . . . , Ci, . . . , Cm I Tf | h, . . . , bi, . . . , b m ) = 5 Clh . . . T c a A . . . 5 Cmbm , (9) 

where T°j, = if ca b (colour-charge matrix in the adjoint representation) if the emitting par- 
ticle « is a gluon and T% = t a a a (colour- charge matrix in the fundamental representa- 
tion) if the emitting particle i is a final-state quark (in the case of a final-state antiquark 
T£g = t®B = — t% a ). Note that the colour-charge operator of an initial-state parton is de- 
fined by crossing symmetry, that is by (Ti)% = t a afi = —t^ a if i is an initial-state quark 
and (Tj)% = t a p a if i is an initial-state antiquark. 

In this notation, each vector \A4 m (pi, ...,p m )) is a colour singlet, so colour conservation 
is simply 

m 

Y J T l \M m ) = U. (10) 

i=l 

The colour-charge algebra for the product (Tj)°(Tj)° = Tj • Tj is: 

Ti • Tj = Tj ■ Ti if i^j; lt = Ci, (11) 

where Cj is the Casimir operator, and we have C, = Ca = N c if i is a gluon and Cj = Cp = 
Tr(N% — 1)/N C = (iVp — 1)/2N C if i is a quark or an antiquark (we are using the customary 
normalization Tr = 1/2). 



3 Singular behaviour at one-loop order 

In this section we recall the known results |T7], [5|, [5| on the infrared-singular behaviour of 
QCD amplitudes at one-loop order. The one-loop coefficient subamplitude M.$(fj?; {p}) 
has double and single poles in 1/e. The coefficients of these poles are universal and can be 
given by the following formula 

|a^V;M)>r S = J (1 W 2 ;M) \mS^{p}))^ + \mS^;{p})) m . (12) 

The contribution A4$ &n on the right-hand side is finite for e — > and, hence, in 
Eq. ( [T2] ) all one-loop singularities are factorized in colour space with respect to the tree- 
level amplitude . The singular dependence is embodied in the factor 1^ that acts as 
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a colour-charge operator onto the colour vector |.M^). Its explicit expression in terms of 
the colour charges of the m partons is 

^W^,.-—^ ^vrW gT.-T, pL_) , ,13) 

where e _ * Aij7r is the unitarity phase (Ajj = +1 if i and j are both incoming or outgoing 
partons and Ay = otherwise) and the singular (for e — ■> 0) function Vj Smg (e) depends only 
on the parton flavour and is given by 



1 1 



VT ns (c) = T?^+7i- • (14) 
e z e 



The flavour coefficients T 2 and 7j are 

1 q — 1 q — , J- g — , 

l q = 1*=\Cf , 7 fl = y C A - jjraty . (15) 

Note that in Eq. (|T2] ) the double poles 1/e 2 are factorized completely and not only in 
colour space. The simplest way to see that is to expand Eq. ([T3|) in powers of e and then 
use the colour conservation relation fllPf) , i.e. J2j^iTj = — Tj. One obtains the result 

/i 2 ; M) = \ E \ 2 E Ti ■ T, + 0(l/c) = - E r f + 0(l/c) , (16) 

that explicitly shows the absence of colour correlations at 0(1/ e 2 ). Nonetheless, single 
poles 1/e are still colour-correlated. 



In the factorization formula ([12]) we have introduced the subscripts rs to explicitly 
recall the RS dependence of the various quantities. Note that the insertion operator 1^ is 
RS-independent. The RS dependence of the one-loop amplitude is completely taken 
into account by that of the tree amplitude JM/$ and by the finite remainder fin . In 
particular, the interference between RS-dependent terms| of (9(e) in and double poles 
1/e 2 in 1^ produces, in general, an RS dependence of M.^ that begins at 0(1/ e). 

The RS-dependent terms of fin are not colour-correlated. Since they are completely 
factorized, we can avoid the colour-space notation and write 

M$*i(A M) = \ (Ett) m^ k ,^ 2 ; M) + f&V; M) + o(e) , (n) 

where the finite coefficients jf a depend only on the flavour of the external partons in Ai m . 
The transition coefficients 7? a that relate the RS mostly used in one-loop computations 
were first calculated in Ref. |[28|| . Using the normalization in Eq. (|T7j ) and assuming CDR 
as reference scheme (i.e. setting ^P DR = 0, by definition), the explicit expressions of 7f a 
can be found in Ref. 



The residual RS dependence on the right-hand side of Eq. fll7D affects only the terms 
of (9(e). In particular, the e-independent function F^\fi 2 ; {p}) is RS-independent [^8], |7 



§Note that the e-dependence of any amplitude M. m is always understood. 
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4 Factorization formula at two-loop order 



The two-loop coefficient subamplitude M$ has 1/e 4 , 1/e 3 , 1/e 2 and 1/e poles. Because of 
the increased degree of singularities, it is not a priori guaranteed that all of them can be 
controlled by a universal factorization formula as in the one-loop case. The main result 
presented in this paper is that such a factorization formula does exist and can be written 
in the following form 

l-M^V;M)> R .a = I {1) (e,^{p})\Mg\^{p})) K& 

+ 4 2 1(^Mp})|A^V;M)>b.s+ l-M^V;M)>aa - (is) 

The factorization structure of Eq. (|i~8"D is not trivial and is somehow analogous to that of 



the collinear limit of one-loop amplitudes considered in Ref. pO 



The main features of Eq. (18) are the following: 



In the first term on the right-hand side of Eq. (0), the one- loop insertion operator 1^ 
of Eq. flT3|) acts onto the one- loop matrix element. The corresponding contribution 
thus contains poles of the type l/e n with n = 1, ... ,4, coming from the single and 
double poles in and ■ 

The second term on the right-hand side of Eq. ( f[q) contains a new colour-charge 
operator 1^ that acts onto the tree-level subamplitude. The two-loop insertion 
operator J*- 2 ** is given as follows 

i ( H(e^ 2 ;{p}) = -i/ (1) (e^ 2 ;M)(/ (1) (e,^;M) + 4vr/3 i 

+ ^[i 1 ~ 2£) (27TA.1 + K) IW (2e, ^ {p}) (19) 
+ H ( H(e^ 2 ;{p}) , 
where the coefficient K is: 

The first and second lineal on the right-hand side of Eq. flT9|) lead to e poles that are 
at most of 4th- and 3rd-order, respectively. Moreover, these two lines control all the 
poles up to 1/e 2 . In fact, the remaining term H^l contains only single poles: 

HW(e,^;{p}) = 0(l/e) . (21) 

The last term, A4$ Rn , on the right-hand side of Eq. (|I~8"| ) is a non-singular remainder 
in the limit e — > 0. 



^Note that on the second line the argument of the operator 1^ is 2e rather than e as on the first line. 
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Note that the two-loop operator /j^ in Eq. (|T9|) depends on the RS only through the 
single-pole contributions in H^. At present, we cannot give an explicit expression for 
H^l that is valid for any QCD amplitude. The particular case of a qq pair is considered 
in Sect. 



Using the factorization formula fll8| ) and Eq. (|19|) , all the coefficients of the poles 
1/e 4 , 1/e 3 , 1/e 2 can be explicitly evaluated in terms of the one-loop operator the first 
coefficient /3q of the beta function and the constant K in Eq. (|20|). 



In particular, the coefficient K allows one to control the double poles 1/e 2 in a universal 
way, that is, independently of the given QCD amplitude. The same coefficient typically 
appears in the resummation program of higher-order logarithmic corrections of the Sudakov 



type ||29|| . This universality |30[ follows from the fact that K measures the renormalization 
(in the MS renormalization scheme) of the intensity of the lowest-order soft-gluon emission. 
The origin of K is thus essentially ultraviolet, which explains its independence of the RS 
of infrared singularities. 

The non-singular remainder _M^ fin is anologous to M^ )Rn in Eq. ([12]). However, its 
RS-dependent part is in general not factorized in colour space and cannot be expressed in 
factorized form as in Eq. (JTT^ 



5 Two-loop amplitudes with two partons 
5.1 Colour structure 

Amplitudes with only two QCD partons have a trivial colour structure to any loop order. 
The colour-space formulae of Sects. || and f| thus become 'true' factorization formulae 
because all the insertion operators are proportional to the identity matrix in colour space. 

In the case of the one-loop insertion operator I^ 1 ', one can use colour conservation as 
in Eq. (|10D and explicitly perform the colour algebra, that is Ti-T 2 = —T\ = —T 2 2 . Then 
one obtains 

T(l - e) \ 2pi -p 2 J 

where the subscript ij denotes the flavour of the two partons: ij = gg or ij = qq, qq, qq for 
the various crossed channels with two fermions. Inserting Eq. (p2|) into Eq. (|l^), one can 
thus rewrite Eq. (|T8|) in an explicitly factorized form. 



5.2 Amplitudes with a qq pair: complete structure of singular 
terms in CDR 

The only QCD amplitude that has been computed so far at two-loop order is the electro- 
magnetic form factor of the quark [J0[ p| . The result of Ref. [|llj and Eq. (|2~2|) can be 
used to check our Eqs. ( |T8| ) and (|T9| ) in this particular case. As a by-product of this check, 
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we can moreover derive the explicit expression of the 0(1/ e) operator ffg for the QCD 
amplitudes with a qq pair or, in general, for those related to them by crossing symmetry. 
We obtain 



1 e -etf(l) 



Ve" Ul27r ' 
, 2pi • P2 



2c 



4e r(l - e) 

1 56 /1fi 

- 7(i) + 3C F K + 5( 2 7c(3oC F - -7i(3 C F - f — - 7C 3 ) CTpCU 



(23) 



where 



17 8c 



7(i) = (-3 + 24( 2 - 48C 3 )C^ + (-y - y C 2 + 24C 3 ) C F C A + (- + — C F T R N f . (24) 

Note that the Cf-part of H^q ie entirely included in the coefficient 7(1). This is exactly the 
same coefficient as controls the virtual contribution (i.e. the term proportional to 5(1 — z)) 
to the NLO Altarelli-Parisi splitting function in the flavour non-singlet sector |3T| . 



We recall that the operator in Eq. fl2"3|) is not only related to the electromagnetic form 
factor of the quark. The same operator applies to the two-loop factorization formula of any 
QCD amplitude with two quarks such as, for instance, the amplitude qq — > 77, which is 
relevant for the hadroproduction of diphotons. Note also that expression (|23|) is valid only 



in CDR because it is obtained by a calculation [] 1| that uses this RS. 



6 Two-loop amplitudes with three partons 



6.1 Colour structure 

In the case of QCD amplitudes with three partons, the one-loop insertion operator 1^ 
is again factorizable in colour space. Using colour conservation (i.e. Eq. flT0D), we have 
2Ti • T 2 = T3 — T\ — T\ and likewise for the other permutations. Then, in the case of 
amplitudes with three gluons, we obtain: 



1 p-^(i) 



ye^y f^e^ 



2pi • P2 



2P2 ■ P3 



2P3 ■ Pi 



(25) 

and, for amplitudes with two fermions of momenta pi,p 2 and a gluon of momentum p$, we 
get: 



1 e -«Ml) f 1 /,;-, 

r VT ( c ) 7T (W F - C A ) ±Z 

) I C F \ 2 Pl ■ p 2 



2 rfi 



(26) 



2p 2 ■ Pz 



2pz ■ Pi 



The fact that the insertion operator 1^ is proportional to the identity matrix in colour 
space does not imply, however, that formula (|18D is exactly factorizable in the three-parton 
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case. In fact, colour correlations can still be present in 1^ ' through the insertion operator 
H^ 2 \ When acting onto the tree-level matrix element , the operator H^ 1 can lead to 
colour transitions between its possible colour states. 

In this respect, qqg and ggg amplitudes behave in a different way. In the qqg case, there 
is only one possible colour state, namely oc t^ a , and thus exactly factorizes. 

In the ggg case, the gluons can be either in the symmetric (Mg gg oc f a b c ) or in the anti- 
symmetric (■Mgl'g oc d a b c ) colour configuration, and the colour structure of has to be 
properly taken into account. 



6.2 qqg amplitudes and e + e — ► 3 jets at NNLO 



Perturbative QCD predictions at NNLO for three-jet observables in e + e~ annihilation are 
strongly demanded by the high experimental accuracy of the LEP and SLC data ||. These 
predictions require the calculation of the corresponding scattering amplitude at two-loop 
order. Since the knowledge of its singularity structure can help to set up and to check the 
calculation, we apply in this section the factorization formula flI5|) to this particular case. 



We are interested in the process e + e~ — > q(pi) + <?(p 2 ) + g{pz)i where Q = pi + p 2 + P3 
denotes the total four- momentum and yij = 2pi ■ Pj/Q 2 are the relevant dimensionless 
invariants. To compute the corresponding cross section at NNLO, one has to evaluate the 



interference \A4 



ggg I (2-ioop) 



between the tree-level and two-loop matrix elements: 



\M 



</ggl(2-ioo P ) = (^fqg I ^Sg) + complex conjugate . 



The analogous interference \M-qq g \ 2 i-\ oop ) at the one-loop level, 



l-M^sla-ioop) = (MfllMgg'g) + complex conjugate , 
was first computed in Refs. [33(1 . 

Using the factorization formulae ( ]T2| ) and ([18]), we can write 



(27) 



(21 



I KA I 2 

K vl <?gsl( 2 -loop) 



2 \M 



(0) |2 

ggg\ 



Re4|-(lml«) ; 



+ \M 



2 Vt» 
99sl(i-loop) rx '^- L qqg 



0(l/e) 



Then, we can use the explicit expression for \-M q q g \ 2 i-i oop \ and Eqs. ( f26|) and (|T 
obtain the final result: 

\Mggg\ (2 

—loop) 



(29) 
and 



2 _ I \yl(0) |2 

ggg Irs i ^4 



+ 



\M 
1 

1 

2! 



2e 3 



[/2(e)/i(e)+67rA,/ 2 (e)-7rA)/ 2 (2e)] 



ff(e) - 7r'fi(e) + 87rA)/i(e) + 12^ + 2L K& f 2 (e) - Kf 2 (2e) 



7^2 /2(e) F{y 12 ,y 23 ,y 13 ) + 0{l/e) 



(30) 



where the function F(y 12 , y 23 , y 13 ) is given in Eq. (2.21) of Ref. |32] and the functions / 2 , f\ 
and the coefficient L Ra are defined as follows 



He) 



-eV(l) 



r(i-c) \Q 2 



(2C F - C A )y£ + C A (y£ + y 23 e 



(31) 



10 




^ 2 



__(2C F + CU) + 8C F - 2y™ - % RS . (33) 



The coefficients 7? a in Eq. ( p3]) parametrize the RS-dependence of the amplitude (see 

Eq. fll?] )). The complete evaluation of the C(l/e)-terms on the right-hand side of Eq. (j30|) 

(2) 

requires the identification of the still unknown operator H q q g - 



7 Summary 

In this paper we have presented a first discussion of the singular behaviour of on-shell 
QCD amplitudes at two-loop order. The complete structure of the infrared singularities 
is described by the colour-space factorization formula given in Sect. The factorization 
formula is universal, i.e. valid for any amplitude, and the singular factors only depend on 
the flavour and momentum of the coloured external legs. At present we can explicitly give 
only the coefficients of the 1/e 4 , 1/e 3 and 1/e 2 poles of the dimensionally-regularized singular 
factors. These coefficients are nonetheless known in a form that is manifestly independent 
of the RS of the infrared singularities. The remaining single-pole contributions, namely the 
operator in Eq. (|19"D, still have to be explicitly evaluated. Owing to their universality, 
they can be extracted from the calculation of few basic two-loop amplitudes, as discussed 
in Sect. |57| for the case of amplitudes with a qq pair. 

Our factorization formula can be useful both to check explicit evaluations of two-loop 
amplitudes and to organize their calculations in terms of divergent, but analytically com- 
putable, parts and finite remainders that can be integrated numerically. In the more general 
context of NNLO calculations of jet observables, our two-loop results can be used to set 
up the integration of tree-level and one-loop amplitudes in such a way as to construct 
process-independent techniques for infrared cancellations. 

Acknowledgements. We would like to thank David Kosower, Lorenzo Magnea and Willy 
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